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The long standing problem of proton-neutron pairing and, in particular, the limitations imposed 
on the solutions by the available symmetries, is revisited. We look for solutions with non-vanishing 
expectation values of the proton A p , the neutron A„ and the isoscalar Ao gaps. For an equal 
number of protons and neutrons we find two solutions with A p = ±A n , respectively. The behavior 
^\ ' and structure of these solutions differ for spin saturated (single /-shell) and spin unsaturared systems 

0\ . (single j-shell). In the former case the BCS results are checked against an exact calculation. 
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CN ! I. INTRODUCTION 



There is presently a revival of the interest on the solutions to the old problem of proton-neutron pairing. This 
' is due to the availability of unstable beams and targets and of large array detectors, which allow the experimental 
study of heavier nuclei with equal number of protons and neutrons. Most of the recent papers include rather complete 
' histories of the subject (see, for instance, refs. (!],§]). 

Specific aspects of the formalism may be found in a review article by Goodman |3|. In particular, the topic of 
the limitations on the form of the solutions which are imposed by the available symmetries, is covered there. Since 
then, most papers on the neutron-proton pairing problem have made use of this important work. However, in the 
present paper we show that the solutions may take more general forms than those found by Goodman. For the sake 
of completeness we reproduce Goodman's arguments in Appendix ^ and indicate the point at which we depart from 
them. 

In section |n| we construct and solve the mean field approximation, starting from a number- and isospin-conserving 
hamiltonian, and applying an appropriate transformation from particles to quasi-particles. As a result of such trans- 
formation the proton A p , the neutron A„ and the isoscalar Ao gaps may take non- vanishing values. We look for 
non-trivial solutions in which the three of them are different from zero. Although the formalism is at least valid for any 
situation involving separable pairing interactions, we restrict the application to nuclei with equal number of protons 
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and neutrons. The nature of the solutions thus obtained is discussed in detail for a single Z-shell (section III) and for 



a single j-shell (section IV). A comparison between exact and approximate results is made in the former case. Since 
significant differences may be found in the behavior of the solutions for spin saturated and unsaturated systems, a 
suggestion is made in the Conclusions (section |v|) about the regions of the nuclear chart in which non-trivial solutions 
may be favored. 

In the main part of the paper we include only general arguments. All the details of the calculations can be found 
in the Appendices [b|, |c] and |5| 

In the present paper we confine the discussion to the problem of the deformed mean field treatment. Improvements 
over this approximation involves the relation between the laboratory system, and the instrinsic system in which the 
deformation takes place. Projection techniques have recently been used to this purpose [0. 



II. THE GENERALIZED BCS TREATMENT 



A. The transformation 



Let us consider a general, number- and isospin-conserving hamiltonian. Assuming that it may be treated within a 
(gauge- and isospin-deformed) field approximation, we define the following transformation to quasi-particles: 



1 




1 np 



(2.1) 



with u 



VW 7 U VW 



real, where v = (proton, ncutron)= (p, n) = (1,-1); A denotes either j or (I, s), while u labels the 



corresponding projections 1 (/i > in (2.1)). The operator c\ Xfi creates a particle in the time-reversed state of c v 



The time-reversed and hermitean-conjugate versions of (2.1) complete an 8 x 
4x4 ones. 

simplifies over the one given in ref. 



transformation, which reduces to two 



The transformation ( 2.1 



on two accounts: 

i) It does not mix operators creating (or annihilating) particles in time-reversed states. As shown by Goodman, this 
pairing only becomes relevant for well shape-deformed nuclei. 

ii) The transformation (2.1) only yields non-vanishing expectation values for the components J = 0,T = l,To = 
v = ±1 and J = l,Jo=T = Oof the tensor {Sj^ To y constructed as a product of two single-particle creation 
operators coupled to spin J and isospin T, with projections Jo and To, respectively. In spite of this limitation, the 
transformation takes into account the full neutron-proton pairing. The case is completely analogous to the familiar 
situation in shape-deformed nuclei, in which the expectation values of the quadrupole tensor are: 



<0o>^0 , < (Q2 + Q-2) >^0, 
<Q±i > =< (Q2-Q-2) >=0. 



(2.2) 



The first line, non-vanishing expectation values are the (lar ge) "order parameters" of the problem, while the ones in 
the second line define the intrinsic system. Equations (2.2) violate the rotational invariance in the intrinsic system 
but not in the laboratory system. We may restore such invariance through the introduction of collective coordinates, 
within an appropriate formalism (see, for instance, Ref. Some of the possible order parameters may vanish if 

there is a subsymmetry still conserved, like for an axially symmetric deformation. As a consequence, the Nilsson 
model takes into account the full, spherically symme tric, quadrupole interaction within the mean field approximation. 

In a completely similar way, the transformation (2.1) may be used to treat neutron-proton pairing problems. In 
particular, we do not need to include any real non-diagonal components in the transformation (2.1), which would yield 
a non vanishing expectatio n va lue for the component T = 1, To = J = 0, since this component plays a similar role as 
those in the second line of (2.2). The fact that in the deformed field treatment of an isovector pairing force there are 
only two order parameters and that, through a convenient orientation of the intrinsic system in gauge- and isospace, 
such parameters may be chosen as the nn and pp gaps (orientation A) is demonstrated in Refs. and 0. Other 
(equivalent) orientations of the intrinsic system are possible such as the one in which the two order parameters are 
represented by the value of the np gap and the equally valued pp and nn gaps (orientation B, cf. Ref. Again, the 
introduction of an additional subsymmetry between protons and neutrons may eliminate one of the two parameters: 
either the difference between the nn and pp gaps in case of orientation A or the nn and pp gaps in orientation B. This 
is consistent with the results obtained in Refs. B and |^] for nuclei with equal number of protons and neutrons, that 
there are two BCS-like solutions with T = 1 pairs: one with no np pairs and the other composed entirely by np pairs. 
Both are degenerate and do not mix with each other. We argue that they correspond to the same state described from 
different intrinsic sytems. Obviously, the physical results obtained are independent of the (unphysical) orientations 
of the intrinsic system, if a proper treatment is applied for relating the values of the magnitudes in the intrinsic and 
laboratory frames of reference]]. 

The orthonormalization conditions are: 



v sv ) = 0. 



(2.3) 



In the present paper we use the transformation (2.1) in order to construct the independent quasi-particle hamil- 
tonian. We also study properties of the vacuum state | >, which is annihilated by the quasi-particle annihilation 
operators a v \^, i.e., 

+ (UnnVpp — u npVpn)Cp Xfi C^ Xfl + (u pp V nn ~ u pnV n p)c < nX ^C < nX - 



1 More details of the notation are given in Appendix 

2 The formalism of Ref. |H) has been recently applied to the isospin case |lC 
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+i(u nn V np M„ p W„„)(c| lAAI Cp A - C pA^ C nA/i) 

+(VppV nn - v np v pn )c\ x ^ nXfi cl Xll cl Xji ) |vacuum > . (2.4) 

In spite of the fact that the calculation is made in the intrinsic system, neither the vacuum energy nor the quasi- 
particle excitation energies are modified by the corrections needed to restore the symmetries (at least to leading order, 
within an expansion on the inverse of the order parameters) . 



B. The solution 



In order to keep the formalism as simple as possible, we consider in this work a separable pairing Hamiltonian of 
the form 



vX fj, 



To 



Jo 



It is useful to parametrize the interaction strengths as 

.9i =.9(1- x) ; .go =g(l + x) ; 



-1< x < 1. 



The matrix elements of the transformation (2.1) are fixed through the following steps: 
i) The diagonalization of the two-quasi-particlc Hamiltonian 



tf (11) = E e -£(4^v) (11) 

v,X \i 



(2.5) 



(2.6) 



(2.7) 



where the supraindex (11) denotes as usual the product of a quasi-particle creation and an annihilation operator. 
Using the shorthand notation S v — SqI and So = Sqq, the gaps are defined as A„ = g\ < S v > and Ao = go < So >■ 
The single-particle energies e v \ — e v \ — X v include a Lagrange multiplier A„ . The diagonalization of the hamiltonian 
(|2.7|) is equivalent to the matrix diagonalization (cf. Appendix 0) 



c„a ^a^Ao A„ 
A p k XfJ .A 
VA:a m A A„ 



-£ p a 






-EnA 









( Upv \ 


) 


\Vnv J 


— EvXfi 


l^nv 1 

Vpy 1 

\Vnv J 



(2. 



that yields the quasi-particle energies E v \^ and the coefficients u vw , v vw as functions of the single-particle energies 

and gap parameters. The factors k Xfl are defined in Table |. 

ii) The fulfillment of the self-consistency conditions [cf. Eq. flB2|)]: 



A^, Q\ ^ ^vw^vw 5 
u>, A,/i>0 



Ao = .9o £ k Xll (u pv v nv + u nv v pv ) 

v,\,[j,>0 



iii) The fulfillment of the number equations [cf. Eq. ( p32| )] : 

w,X,fi 



The vacuum energy is given by the expectation value (< H > 

w, 



A 2 4- A 2 

" ~2^ e ^2^ U ™ g(l _ X ) 



A 2 

^0 



3(1 + x) 



(2.9) 



(2.10) 



(2.11) 



In this paper we are interested in situations with the same number of protons and neutrons. If we further assume the 
same single-particle spectrum for both kind of particles (e p \ — e n \), the symmetry of the problem requires A 2 = A 2 
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but leaves open the relative sign of A p and A n . Therefore we expect two non-trivial solutions with A p = ±A n ^ 0, 
A 7^ and two trivial solutions^ with A p — A n = 0, A ^ and A p = A„ ^ 0, A = 0. 

To explore the characteristics of these solutions, in what follows we work in a model space consisting of one 
single shell. In this case, the single-particle energies e v \ disappear from the problem and all energies become simply 
proportional to the strength parameter^]. 

We have the option to work with an Z-shell or a j-shell. There are at least two advantages for considering the case 
of a I- over a j-shell: i) The pair of particles may couple to orbital angular momentum L — 0, even for J = 1. ii) 
The pairing problem is amenable to a numerical solution and we may check the BCS results against the exact ones. 
Nevertheless, we also carry the calculation for the case of a single j-shell, in order to study the effect of the pairing 
hamiltonian in a spin non-saturated system. 

III. THE CASE OF A SINGLE i-SHELL 
A. The BCS solution 

The solution with A„ = — A p is the one allowed by Good man but for the fact that the transformation to quasi- 
particles displays non-diagonal coefficients u vw ^ [Eq. (|Cl[)]. However, since the eigenvalues of the Hamiltonian are 
degenerate, it is always possible to make linear combinations of the solutions corresponding t o th e quasi-protons and 
quasi-neutrons (i.e., to the two values of v) and impose the condition u np = u pn = [cf. Eq. (|2.g|)]. 

However there is a sharp limitation in the domain of validity of this solution given by the self-consistency condition 
) , which implies x = [(cf. Eq. ( |2.6| ))]. Moreover, the allowed interval around this value has zero width. 



[Eq. (C2) 



Figure llrepresents the vacuum energies as a function of x for two different values of the occupation parameter rj 
(see TableH). The energy W gs associated with the solution having A p — — A„ lies exactly at the crossin g po int of 



the energies corresponding to the two trivial solutions (Aq = or A p = 0). Since the solution given in (CI) only 



determines the total gap y A p + ^Aq, but not the ratio between the isovector and isoscalar gaps, the existence of this 

solution may be interpreted (in the present case) as a manifestation of the degeneracy associated with the crossing, 
without further physical meaning. 

The two trivial solutions exist for any value of the parameter x. 



The solution with A„ = A p requires the numerical treatment of two simultaneous equations [i.e., ( |C7|) and (C8)], 
depending on the two variables 

C = e/A p and 7 = A /A p . (3.1) 



The t ask is simplified by the existence of transformations that leave invariant the system of equations [Eqs. ( C10| ) and 



flCllD ]. Thus it is sufficient to solve the system in one quadrant of the cartesian system determined by the strength 
parameter x and the occupation parameter 77. 

In Fig. |l|(a) there appears an upper value of |x| for which there is a solution. In Fig. [|(b) there is also a lower value 
of I x I . In the following we explain this behavior. 

According to Fig. El the vacuum energy W gs of the solution with A p = A„ lies always higher than the energies 
of the trivial solutions and approaches these preceding values at the extremes of the domain in which it exists. For 
instance, for negative values of x, the energy of the solution with A p = A n becomes the energy of the highest trivial 
solution, namely the one with a non-vanishing isoscalar gap. This fact suggests that one limit of the domain in the 
(x,r/) plane is obta ined from solutions such that 7—5-00 and £ — > 00. In such a limiting case the following curve can 
be found [cf. Eqs. dcl^: 



/l + si (-l<:n<0). (3.2) 
V 1 -Xl 

The procedure also yields the (vanishing) value of A p ; the (non- vanishing) value of Ao; and the vacuum energy W gs 
on this curve. We may verify that they are equal to the results of the trivial solution with A p = that are given in 
Table || 
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The relative sign of A p and A n is immaterial for the two trivial solutions. 

1 

2 



4 We take g — | in the numerical calculations below 
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The limiting curve for positive values of x\ is easily obtained from (3.2) by means of the transformation ( |C10| ) 
(Note that in this case the gap A p -> A /V2 and A -> [cf. Eq. flCl^)].) 



m = 



1 — X\ 
1 + 



(0 < xi < 1). 



(3.3) 



The curve (3.3) that limits the allowed region is shown in Fig. g. 

According to Fig. [I], the solution with A p = A„ does exist for rj = —0.80 in an interval around x — 0, but it does 
not exist for 77 = — .46 in an interval enclosing the same point. In order to explain this behavior we construct a second 
boundary for the domain of validity by studying this so lutio n in the limit £ — ► and 7 — ► \2. The locus at which 
this condition is satisfied in the (x,rj) plane is [cf. Eqs. ( p!7 )]: 



(3.4) 



Equation (3.4) indicates that the curve does not exist for values of the occupation parameter 77 > i. Therefore, for 



such values of 77, the region of validity is only limited by Eqs. (3.2) and (3.3). For values r\ < h, the allowed solutions 



lie in the interval X2 < x < x%. This interval gets narrower as we approach the half-filled shell and, in fact, it vanishes 
in the limits x\ — > x 2 — > 1 (Fig. ||). These predictions are consistent with the numerical results presented in Fig. |l|. 

The absence of solutions within the interval < x < X2 may be interpreted as due to the presence of unstabilities 
in the generalized BCS solution: the quasi- neutron energy E n vanishes along the curve (3.4). 



B. The comparison with the exact results 



The Hamiltonian of Eq.(2.5) is invariant under the transformations of the group SO(8). Thus it is solvable by an 
appropriate choice of the infinitesimal operators pi| , |i2| . The eigenstates \k > of the Hamiltonian can be expressed in 
terms of a basis labeled by a chain of irreducible representations of 50(8) and its subgroups, namely, 



\k;SM s ,TM T >=Y J c n ) (SM s TM T )\Q,X, [n,n,0},SM s ,TM T >, 



(3.5) 



where 51 = 21 + 1 labels the representations of 50(8) for the case of shell seniority zero, A = f2 — {A/2) (A = N p + N n ) 
is related to the representations of SO(7), while [n,n,0] labels the representations of the group 5t/(4) which is 
homomorphic to 50(6). The group 5?7(4) can be readily decomposed in SU(2) x SU(2) so that the quantum 
numbers 5, Ms, T and Mr are sufficient to completely specify the states of 50(6). 



We have obtained the solutions of the Hamiltonian as a function of the parameter x introduced in (2.6). The 
corresponding spectra for 5 = and T = are shown in Fig. ||. The value I = 20 has been used. It is sufficiently 
large to allow for a clear manifestation of collective effects. 

If A = 4 << Q, (77 w —1), the exact solution displays only two eigenvalues. The two trivial solutions, Ao = 
or A p = 0, correspond to the lowest eigenvalue for x smaller and larger than zero, respectively [FigJ|(a)]. The 
approximate solution is somewhat less bound than the exact one, as to be expected. The correspondence is inverted 
for the highest eigenvalue. However there is a small interval around x ~ in which the exact eigenvalues display some 
curvature. While the highest eigenvalue may be correlated with the solution with A p = A„, no such correlation can 
be made for the solution A p = — A„ in the case of the lowest state. 

For a larger number of particles (A = 16, 77 = —0.80), there appear more eigenvalues [Fig. [|(b)]. However, as in 
Fig. ||(a), there is a correspondence between the trivial solutions and the lowest and highest eigenvalues, at least 
for values of \x\ > 0.2. Moreover, for this case, a strong correlation can be made between the eigenstate with the 
highest eigenvalue and the solution with A p = A„ throughout the interval —0.2 < x < 0.2. On the contrary, the 
lowest eigenvalue continues to be represented by two straight lines with a relatively sharp crossing at x « 0, which is 
consistent with the absence of a solution with A„ = — A n . 



IV. THE CASE OF A SINGLE j-SHELL 



In order to solve the self-consistency and number equations [(D3), (|D5|) and (D16)] we must select a value for the 



quantum number j. The value j 



21 



is large enough to allow for a clear display of collective effects. 
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The vacuum energy W gs of the two non-trivial and of the two trivial solutions are represented as a function of the 
strength parameter x in Fig. ||. The occupation parameters 77 (see Table [|) are taken to be ry=-4/ll and in Figs. || 
(a) and (b), respectively. 

At variance with the case of a single /-shell, the solution with A p — — A„ exists within a finite domain for comparable 
isovector and isoscalar strengths (x rs 0). Its vacuum energy becomes the lowest energy within this domain. At the 
extremes of the allowed interval, Goodman's solution merges with the most favored trivial solution corresponding to 
the same value of the parameter x. This behavior is common to the two values of the occupation parameter. 

The non-trivial solution with A p — A n also exists for the two occupation parameters. As similarly to the case of 
a /-shell, this solution always has the highest energy. The x domain of existence is somewhat broader than the one 
associated with the solution for A p = — A n . At each extreme of the allowed interval the solution merges with the 
most unfavored trivial one corresponding to the same value of x. 

According to Fig. [|, the vacuum energy of the solution with A p = A„ is a continuous function of x for 77 = 4/11 
and a discontinuous one for 77 = (half-filled shell). In the second case the energy only exists for certain intervals of 
x and is a linear function of x within each interval. 

The case 7; = is treated analytically in Appendix D 2 a . Since the Lagrange multipliers vanishes, the number 
equation is trivially fulfilled. The discontinuities may be tra ced back to the vanishing of the quasi-neutron energies 
E n for certain values of the asymmetry parameter 7 (Eq. (DC)). Such discontinuities determine closed intervals within 
which the self-consistent conditions yield the parameter x as a conti nuo us function of 7. Within each interval, 7 lies 
between the roots corresponding to two consecutive values of m [Eq. (D9)] . The intervals may be labeled by the value 
of fc = m + i, where m is the smallest magnetic quantum number. There are j — \ closed intervals. The f uncti on ^(7) 
may be inverted and thus an analytical function 7(2;) is obtained within each interval [Fig. || and Eq. (D12)]. Note 
that there are several values of 7 which are compatible with a single value of x < . The opposite situation appears 
for x > 0, where there may not be any value of 7 associated with certain values of x. 

The vacuum energies are given in (Eq. D13). Within each interval k, they confirm the linearity with x that 
was found in the numerical solution displayed in Fig. |[ The slope of the lines W gs (x) increases with the value of 
k = m + \. In the first interval the solution is degenerate with the trivial solution Ao ^ and extends within the 
interval — | < x < (to leading order in 1/j). 

We also obtain an expression for the extremes x± limiting each interval (Eq. D14). 

In addition to the j — h closed intervals, there are also open intervals for small and large values of 7. For 7 < 
1/ y/2j(j + 1) and 7 > ^ j/[2(j + 1)] the solutions collapse to the points x = =Fl, respectively, and the energy vanishes. 



V. CONCLUSIONS 



We have studied a generalization of the BCS approximation that takes into account both the isoscalar and isovector 
pairing. Although for separable forces the solution is always amenable to the diagonalization of a 4 x 4 matrix 
equation, in this paper we have only considered the case of identical number of protons and neutrons and identical 
single-particle energy levels for both kind of nucleons. We perform the calculation in an intrinsic system in which the 
non- vanishing gaps are the proton gap A p , the neutron gap A„ and the isoscalar gap A . We label the solution as 
non-trivial if none of the three gaps vanishes and as trivial if either the isoscalar gap or the neutron and proton gaps 
equal zero. We expect two non-trivial solutions with A p = ±A ra and two trivial solutions. The non-trivial solution 
with the minus sign corresponds to the one allowed by Goodman's argumentation. 

We have studied the vacuum energy of the four solutions in the case of a single I- and a single j-shell. We found 
that the two non-trivial solutions behave differently in spin-saturated than in spin-non-saturated systems. 

In the case of a /-shell, the solution with A p = — A„ does not exist. More precisely, it only exist for the value 
x = of the strength parameters, within an interval of zero width (Fig. |l|). The vacuum energies of the two trivial 
solutions also cross at this point and have the same value as the non-trivial one. The exact calculation confirms the 
non-existence of a solution with A p = — A n , since the lowest energy is well represented by two straight lines crossing 
rather sharply at x = 0. 

On the contrary, the solution with A p — —A n exists within a finite interval in the case of a j-shell. In this interval, 
starting at x > —1/7, it becomes the lowest state of the system. As x tends to —1/7 from above, the non-trivial 
solution tends to the trivial one with Ao = 0, which becomes the lowest state in the interval — 1 < x < —1/7. The 
non-trivial solution tends to the trivial solution with an isoscalar gap at the other extreme of its region of existence. 

The solution with A p = A n exists and has similar properties for both kind of shells: it always displays the highest 
vacuum energy within its regions of existence. For sufficiently large values of |a;| it merges with the highest trivial 
solution. In the case of an /-shell, we obtain expressions for its domain of existence in the plane determined by the 
strength x and the occupation number 77. The energy of the quasi-neutrons vanishes and the system becomes unstable 



G 



along one of the two curves limiting this domain. The pattern is confirmed by checking the BCS results against exact 
calculations: it is possible to identify the highest exact eigenvalue with this non-trivial solution. This general picture 
is very similar for the case of a j-shell. The main difference appears for a half-occupied shell (r] = 0), where the 
regions of validity consist of successive intervals which properties may be predicted analytically. 

Concerning the observation of the effects discussed in the text we may conclude that the experimental search for 
(large) non-trivial solutions should focus on nuclei such that both valence nucleons mainly occupy the j-level which 
becomes integrated to the lower harmonic oscillator shell. 
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APPENDIX A: SOME SYMMETRY CONSIDERATIONS 



In order to reach the point where our results depart from Goodman's, in this Appendix we fol low its presentation 
of the consequences of the symmetries involved H. In Goodman's notation, the transformation ( |2.l[ ) reads: 



U -V 
-V* u* 



(Al) 



where (c*) = 



,t 



'nXfi ' pA/i ' nXfi 



and 





Ilpp 


lU n p 










lUpn 


















Upp 










lUpn 





u 

The generalized density matrix is defined as 



V 









Vpp 


IV n p 








IV pn 






IV n p 








ZVp/}, 











(A2) 



P 

it 1 



i 



where transpose matrices are indicated by a tilde and the p and i matrices satisfy the relations 

p = p f ; i=-t. 

The requirement TZ = TZ 2 is equivalent to: 



P-P 



tt" ; pt — t p. 



(A3) 
(A4) 
(A5) 



Time reversal con servation implies the existence of relations between the matrix elements of p and t. If in addition 
the restrictions (A4) are applied, one obtains p and i matrices of the form 



/ Ppp Pp 

P- 



\ 

Pnn 

Ppp p pn 

\ Ppn p nn J 



pn rnn 




t = 



( o 



t 



l'l> 



tpp tpn 



pn tnn 



V 



-t*pn 
-inn j 



(A6) 



with the diagonal terms p pp , p nn , t pp , and t nn being real numbers. Moreover, the form of the transformation (2.1) has 
the additional consequenc e th at p pn and t pn are purely imaginary. 
Now the requirements (A5) imply 



rpp 
2 

Pnn 



I 2 

on I j 



I Ppn I ipp "H |ipr, 

Pnn — Pnn — |Ppn| — i n n "I" |ipn| ; 
Ppn (1 ipp inn) — ipn (inn ~t~ ipp): 
Ppn (ipp inn) — ipn (Ppp Pnn) • 

Let us consider the expectation value of the rising isospin operator 

-i > = --7= (^"p + ^«p) 



(A7) 



< T 4 



jm>0 



= --L £ (Pnp 



+ P*np) 



(A8) 



jm>0 



Thus this expectation val ue is always zero for transformations which yield purely imaginary values for the non-diagonal 
matrix elements p np (as ( |2.l| ) does). 

If N n = N p and p np — 0, Eqs. ( |A7[ ) indicate that i pp = —inn (Goodman's solution). However, although such 
solution may exist, there might be also another solution with p np ^ and t pp — t nn . 
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APPENDIX B: THE CONSTRUCTION OF THE GENERALIZED BCS BASIS 



In the following we define the pair operators S v , So in the A,/i notation. 



A,ai>0 

Si = -i E M<&^U - 4^). (Bl) 

A,/i>0 

The equivalence with in the Is and jj coupling schemes is given in Table |. as well as the coupling factor k\^. The 
pairing operators ( Bl ) and the sin gle particle operators may be expressed in the space of the quasi-particles using the 
inverse of the transformation (2.1). The result for the vacuum terms is: 

*C S v ^> — ^ u vw v vw , 

A,/x>0 

< Sq > = ^ k\^(up V v nv + u nv v pv ) , 



while the two quasi-particle components appearing in the Hamiltonian ( p.7[ ) yield 

+i 22 { u vpVvn+u vn v vp ){a px ^a n \^ — a pXji a n \p) +h.c., 

\,fi>0 

(Sq + So)^ — -2 E k\^(u pv v nv + u nv v pv )al Xfi a v x^ 

^ ^ k\^(UppV nn -f~ U nn Vpp -\- Up n V n p -\- UmpVpyi) 
\,ft>Q 

x Ka m "™a m - al Xfl a nXfl ) + h.c. (B3) 
Finally, the components of the single-particle Hamiltonian are written 

\ (ii) 

/] e v\ 22 ^AjU^V I = ^2 6vX /2( U vw ~ V vw) a wXtJ, a w\ti 

am ) v,\ w,fj, 

v.X fi>0 



X OpA^V - a pA/i a nAp) + h.C. (B4) 



Using (B3-B4), the diagonalization of Eq. ( |2.7| ) yields the two equations: 

w w 
-\-2/\ok\fi(Up V V nv + U nv Vp V ), 



-\- Aok\p^ (llppVnn ~t~ ^nn^pp ~t~ ^np^pn ~t~ ^pn^np) i (^^) 



which are equivalent to the matrix equation ( |2.8| 
The (positive) eigenvalues are: 
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(B6) 



D = 



^(e 2 pX + - el x Air + 4Ag^ ((e pA - £ „ A ) 2 + (A p + A„)^ 



where we may (arbitrarily) assign the highest energy (v—1) to the quasi-proton. The eigenvalues depend on the 
magnetic projection fi through the factor k\ a - Therefore, according to Table |[ they do so in the case of jj coupling, 
but not for Is coupling. 

In this paper we are interested in the situation e n \ — e p \ {— e\) and N p = N n The expression for the eigenvalues 
simplifies to: 



E, 



vXfj- 



e 2 x + A 2 p + A 2 k 2 Xfi + vA Q k Xfl (A p + A n ). 



(B7) 



The eigenvectors of (2.8) yield the coefficients of the quasi-particle transformation (2.1). If the conditions leading 
to (B7) are valid, we obtain the (unnormalizcd) coefficients 



-A k Xfi (A P + A n + v^A 2 k 2 Xu + {A p - A n ) 2 

+ (A n - A p + v^4A 2 k 2 Xn + (A p - A r , 



Vpv — A k\ u (E v \ p — e\), 
1 



(E vXn - e x ) [An - A p + vJAA 2 k 2 Xu + {A p - A„) 



(B8) 



APPENDIX C: THE SOLUTIONS FOR THE CASE OF A SINGLE /-SHELL 



1. Solution with A„ 



-An 



According to Eqs. (|B7|) and (|B8|) 

E V = E= \/e 2 + A 2 



<• I A + vA p I E + i 



E \E-e 
1 lA + vAp (E-e 



I A- vAn E 



E 



E 



E 



A - vA, 



E 



The self-consistent equations (2.9) and the number equations (2.10) yield 



The vacuum energy (2.11) is 



9i= 9o = 9 = 



2E 
21 + 1 



V = -E- 



E-e 
E-e 



E 



(CI) 



(C2) 



W, 



3 s 



,(2Z + l) a (l-^). 



(C3) 



2. Solution with A p = A r , 



E v = ^e 2 + {A p + vA /V2) 2 = A p yJ ( 2 + (1 + V1 
_A p + wA /\/2^ _ vA p + A /V2^ 



2y/E v (E v -e) ' 2 x /E v (E v -e) : 

1 lE v -e v \E v -e ^ 



Z V Xv-u Z V -&v 
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where £,7 are given in (3.1). The self-consistency equations (2.9) and number equations ( 2.1 0| ) read 



1 9l ^r n { Ep + En + ^k {En Ep] 



9o 



21 + 1 

AEpE n 



E p + E n 



V2 



(E n - E p ) , 



E p + E n 

V = — e - 

2E p E n 



(C5) 
(C6) 
(C7) 



A combination of the the Eqs. JCq) and (|C6|) yields an expression for the parameter x measuring the relative strengths 
of the two pairing interactions fl2.6|). 



(E n -E p )(^ 



2 VI) 

7 ' 



V2\ ' 



2(Ep + E n ) + (E n -Ep)(^ + ^) 



(C8) 



In principle we should simultaneously solve th e th ree eq uati ons (C5)-(C7) in order to obtain A p ,A ,e as functions 
of gi,go,rj. However it is easy to verify that (C£) and (C7) depend only on the two variables 7 and £. Therefore 
it is simpler to solve first the sys tem of t wo e quations and subsequently obtain the values of A p , Ao through an 
independent combination of Eqs. ( |C5|) and (CE). This last step yields 



A p = ^^±Hx 

AEpE n 



\ 



(1- ){E p + E n + ^={E n - E p )j lE p + E n + ^(E n - E p )Y 



Ao = A p 7. 



These gaps may be used to obtain the vacuum energy W gs (cf. ( |2.1l|) ). 



(C9) 



a. Symmetries and domain of existence 



The system of Eqs. (JCVD and (C8) remains invariant under the two following transformations 

CV2 



7 



c 



-x; T}-> T), 



7-^7; c -> -C ; x -> x ; v^-v- 



As a consequence of the invariance under the transformation (C1C) 



A p -» -| ; A -> V2A P . 

Thus this invariance implies 

W gs (x,rj) -> W ga (-x,ri). 
Since A p and Ao remain invariant under the transformation ( |Cll| ), 

W gs (x,r{) -> W gs (x,-rj). 



(CIO) 
(Cll) 

(C12) 

(C13) 



(C14) 



In order to find the region of allowed solutions in the (x, rf) plane we discuss the two limits Q — > 00 and £ — ► 0. 
For the first limit, Eq. (C7) requires that 7 = (3^^ (with /3 DO =constant). We obtain 



lim x 

£ — >oc 



ft 



2 + ft 



lim 77 



(C15) 
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The consistency between these two equations yields the curve r)i(xi) (^^) limiting the region of allowed solutions. 
Let us calculate now the gaps and the vacuum energy on this curve. From Eqs. ( |C9| ) we obtain 

lim A p = 0, 



lim A = (21 + 1) X J^~x\ ; (-1 < X! < 0) 

C^OO y/2 V 1 + Pic 

= (21 + 1)X (l+xi),Jl^$, (Gig) 



which implies 



W ga = -(2Z + 1) ^ 2 " (l-r^f) 

= -(2^ + l) 2 ^(l-^). (C17) 

In order to obtain the energy for positive values of x\ we apply again the transformation ( |C10| ). According to ( C13 ), 

= -(2« + l) 2 ^(l-ry?). (CIS) 



Let us consider now the limit £ — > 0. Similarly to the previous limiting procedure, we require 7 — > \/2/?oCi with /3o 
constant. The limiting procedure yields 

lim a; = — ^° ; lim 77 = ■ (C19) 

The last equation implies that (3q is a real number in the interval — \ < r\ < i. Within this interval we obtain the 



curve (3.4), which is complementary to fl3.2j) and (3.3) 



3. The two trivial solutions 

The results are summarized in Tabic |l[ 

APPENDIX D: THE SOLUTIONS FOR THE CASE OF A SINGLE j-SHELL 

The pairing operators used in Section III and in Appendix ^ are given explicitly in the second column of Table |. 
They may also be written in the jj coupling scheme, 



,001 u 21 + 1 r t ti£=o.J=o.r=i _ V2j + 1 r t hJ=0,T=l 



K 10 -iJ 2l + 1 r r 1yti£=o.J=i.r=o 



. (2i+ + l)(i+ + l). 



^ irH C' , 1 C 



t J lJ=l 



Gj 



/ (2j- + l)j- t t J=1 
V 6(^ + 1) LC w- C "J- Jm 

J(^ + + i)(2r + i) / [e t t, ]r _ [c 



12 



where j = I ± \. In this Appendix we consider the case of a single j-shell which may be derived from the previous 
Z-case by including a large spin-orbit term in the Hamiltonian. In such case we may keep only the contribution 
to (SjJ M y corresponding to the valence shell. We multiply this contribution to (S^/o)^ by a factor (chosen for 
convenience) that yields the expression given in the third column of Table |. 



1. Solution with A p = — A„ 



According to Eqs. (B7) and (B8) 



p i(i + i)' 



v A + vA p ( E + e 



E 



E-el ' nv ~ 2 



1 lA-vA p (E + e 



E 



E-e 



1 A + vA„ (E-e 



2 V E 

The self-consistent conditions and number equation are: 



E + e ' 2V E \E + e 



(D2) 



es 



i = \gis° = \90s2 ; v = , 2j + , ■ 



E 



Si = > i III , 

m>0 



/(./ - d; v^p + ( l) E n 



(D3) 



In the limit Aq — > 0, the self consistent conditions yield the limiting value x > —1/7. 



2. Solution with A„ = A„ 



E 



v — \ 



e 2 + A p + vA a m,l- 



A p + vA a m^-^ 

2y/ E v (E v -e) 

1 lE„-€ 



3{3 + l)) 1 

vA p + A m ' 2 



2^/E v (E v -e) ' 



2 V JS« 



2 V E v 



(D4) 



The self-consistency equations (2.9) and number equations (2.10) read 



es 



1 = (so + 7*0 = ^.9o(^ + a 2 ) ; »? - 2 . + , • 



(D5) 



where the Si are given in Eq. (D3) 



The analytical treatment of the solution in the middle of the shell (r) = e = 0) 



According to Eq. (D4), the quasi-particlc energies arc (m > 0): 



E p — A p 1 + 7TOi/ — 



1 — jrriA — 



3(3 + 1) 



(D6) 
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which yields the result s of Table HI. With these results we may evaluate the expressions appearing in the self- 
consistency conditions (D5) 



s + jsi = 



p m>0 



Si + 7S2 



2 

a7 



E 

m>0 



(1 - y m ), 



(D7) 



where the values of y m are also given in Table III. Therefore, the relative strength x has the value: 

-V + if + E m >0 Vm (7 + 2-) 



(D8) 



(j + !) 2 +E m>0 S/m(jVT-2m 
which clearly displays the discontinuities to be expected due to the form of E n at 



1 



iO' + i)' 



(D9) 



In the interval k (k 



|) defined as 



2k- 1 



> 7 > 



y/2jQ- + 1) 
2k + 1 



Eq. (D8) is equivalent to 



7Vi(j + 1)/2 - ((i + ^) 2 - fc 2 ) 



7VJ(J + l)/2 + ((J + 5) 2 - fc2 ) ' 
which may be inverted so as to yield the ratio 7 as a function of the relative coupling strength x 

(j + i) 2 -fc 2 



x) k^j(j + l)/2 



(D10) 



(Dll) 



(D12) 



The value of 7 is represented in Fig. || as a function of x. The discontinuities are apparent. Note also that for x < 
there are several values of 7 which are compatible with a single value of x. The opposite situation appears for x > 0, 
in which case there may not be values of 7 associated with some values of x. 

Within the same interval we may express the vacuum energy W gs as a function of x 



W gs =-»-l{l-x)k' -(l + x) 



(D13) 



which displays the linearity in x within each interval. This behavior was already found in Fig. |] and it is to be 
expected for degenerate shells. 

The initial and final values x±, for each interval k, may also be given as a function of k 



x± 



fcj(j + l)-(2fcTl)(Q- + ^ 2 -fc 2 ) 
kj(j + l) + (2kTl)((j + k) 2 -k 2 ) 



(D14) 



Let us study now the behavior: i) before the first interval starts at m = | and ii) after the last one ends at 
m = j. Note that 7 not only decreases from one interval to the next one (D9) but also within each interval ( D1C| ). 
The corresponding values are given in Table IV. The value x = —1 (+1) becomes associated with a predominant 
isoscalar (isovector) pairing correlation. This is consistent with the fact that the solution with A p — A„ merges, at 
the extremes of the s-interval, with the most unfavored trivial solution having the same value of x. 
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3. The trivial solutions 

The solution with A = is the familiar one with 



E= ^/e 2 + A2; e = - V E; 



W gs = -^(j + \) 2 (l-v 2 )- (D15) 
The case with A p = requires to solve the equation 

v = — K-trY^ 1 (Die) 

With the value of the ratio e/Ao thus determined, one obtains the gap and the vacuum energy. 

A ^-^-V - ^ ; W gs = -^. (D17) 

The solution becomes simplified for the case 77 = e = 0, namely, 



50 , 1,2 / 2 



A °( e = °) = T W + a) V, + i' 

^(-0) = -f^- (D18) 
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FIG. 1. The vacuum energies W gs for a single /-shell as a function of x, (a) for 77 = —0.80 and (b) r\ = —0.46. Dotted lines: 
trivial solutions. Full: solution with A p = A„. Big dot: solution with A n = — A p . 
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FIG. 3. The energy spectrum as a function of x for f2 = 41. (a) A = 4 and (b) A = 16. Dot-dashed lines: exact solutions. 
Dotted: trivial solutions. Full: solution with A n = A p . 
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A 


Z 




At 


m, cr = ±i 


m 


At > 


— Z < m < I, a — | 


m > 




1 


™ / 2 


si 
ft 

^0 


-E,V^ r [^ t K JW1 




t? 


w P 1 

2i + l 1 


N P _ 1 

3+¥ 



TABLE I. The equivalence of the notation A, At in the Is and jj coupling schemes, the definition 
of the pairing operators Sl and of the occupation number r\ (—1 < r\ < 1). 



A p 


I fll (2Z + l) v /l-»? 2 





A 





^ So (2Z + l) v /l-'? 2 


e 




-77.E 


E = E p = _B n 


iffi(2Z + l) 
-i 5l (2Z + l) 2 (l-r, 2 ) 


§So(2Z + l) 
-i 5 o(2Z + l) 2 (l- ?7 2 ) 


w gs 



TABLE II. The solutions with A = (second column) 
and A p = (third column) for a single Z-shell. 





i > my/ jU+1) 


1 < 7™^ . ( 


— + — 

E p ^ E n 


2/A p( l--V) 


-W^ /Ap(1 -^ } 


_L J_ 

Ep E n 




2/A p (l-^) 


Um 


1 






TABLE III. The values of -J- ± -gj- and of j/ m as functions of 7 



7 






E m >n m ^ 




y/2j(j + 1) < 7 < OO 











-1 


o<7< y^- 


1 


j + h 


10' + 1) 2 


1 



TABLE IV. The behavior of the solution for very large and 
very small values of 7 
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